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ABSTRACT 

For  systems  that  are  open-loop  stable,  there  is  a  class  of  feedback  controllers 
that  have  the  property  that  the  closed  loop  system  is  stable  and  remains  stable 
in  case  actuator  outages  occur.  Properties  of  a  special  subclass  of  these 
controllers  are  discussed,  .i 
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jj  1  concerning  the  Integrity  property  of  controllers  that  have  gains  which  are 

_A  determined  from  Liapunov  rather  than  Riccati  eouatlons. 
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INTRODUCTION;  Consider  the  linear  controllable  system 
x  -  Ax  +  Bu 


where  A  Is  stable,  l.e.,  the  eigenvalues  of  A  have  negative  real  parts.  The 
class  of  state  feedback  regulators  of  the  form 

u  =  -BTPx  (2) 


where  P  satisfies  the  Liapunov  equation 


PA  +  ATP  +  Q  «  0,  Q  >  0  (3) 

L 

with  (A,  Q9)  observable  are  of  special  Interest  because  the  closed  loop 
systems  with  related  regulators  of  the  form 


u  *  -LBTPx,  L  *  ® 


(4) 


are  stable.  Such  regulators  may  be  considered  to  possess  Integrity  with 
respect  to  loss  of  Imputs,  that  Is,  stability  of  the  closed  loop  system  Is 
maintained  when  one  or  more  Inputs  Is  set  to  zero.  This  situation  of  loss  of 
Inputs  can  be  represented  by  an  admissible  L  In  (4).  For  example,  the  loss 
of  the  first  Input  may  be  represented  by  taking  L  to  be 


The  class  of  regulators  defined  by  (2)  and  (3)  Is  a  subset  of  the  class  of 
optimal  state  feedback  controllers.  This  fact  and  the  proof  of  the  stability 
for  regulators  of  the  form  (4)  are  given  below,  followed  by  other  results  aimed 
at  characterizing  the  class  of  regulators  of  Interest. 
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Optimality  of  Regulators  given  by  (2),  (3)  and  (4):  The  regulator  given 

by  (2)  and  (3)  Is  the  optimal  regulator  for  (1)  with  respect  to  the  performance 

Index 


00 

J  *  J (xT(Q  +  PBBTP)x  +  uTu}dt 


for  P  satisfying  (3)  and  A  being  stable,  P  >  0,  so  that  Q  +  PBB  P  >  0. 

The  optimal  regulator  for  (1)  with  respect  to  (6)  Is  given  by 

u  *  -BTPx  (7) 

A 

with  P  being  the  positive  definite  symmetric  matrix  solution  to  the  Rlccatl 
equation 

PA  +  ATP  +  Q  +  PBBTP  -  PBBTP  «  0  (8) 

Thus,  P  *  P,  and  the  controls  given  by  Cl)  and  (?)  are  the  same.  Now  conifer  the  control 
given  by  (3)  and  (4).  If  L  >  0,  this  control  Is  optimal  with  respect  to  the 
performance  index . 

J  *  j{xT(<?  +  PBLBTP)  x  +  uV’uJdt  (9) 


x  +  uTL“^u)dt 


since  the  optimal  control  Is  given  by 

u  *  -(L"1)"1  BTPx  =  -LBTPx 

where  0  *  PA  +  ATP  +  (Q  +  PBLBTP)  -  PBU"1)"1  BTP 
*  PA  +  A^P  ♦  Q 


If  L  Is  singular,  the  control  v  *  Tu  Is  optimal  with  respect  to 


J  -J{xT( Q  +  PBLBTP)  x  +  vTv}dt 


‘ 1  **  c  t*  i.  C*J  t>' 0  T* 

u,  :  u  * 

^  -  i  *4k  t  J  {  *' 


:  r ' 
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for  the  system 


X  »  Ax  +  BLTTV 

(13) 

where  TTT  *  L+  (the  pseudo- Inverse  of  L) 

(14) 

since  v  *  -  (BLTT)^  Px 

(15) 

and  0  *  PA  +  ATP  +  (Q  +  PBLBTP)  -  P(BLTT)(BLTT)TP 

=  PA  +  ATp  +  Q  +  pblbtp  -  PBLTTTLBTP 

=  PA  +  ATP  +  Q 

(16) 

Thus,  the  closed  loop  system  given  by  (13)  and  (15)  Is  stable,  l.e. 

x  *  (A  -  BLTT(BLTT)TP)x 

*  (A  -  BLTTTLBTP)x 

*  [A  +  B(-LBTP]x  (17) 

Is  stable.  But  (17)  Is  the  same  closed  loop  system  as  that  obtained  using  the 
control  (4)  In  the  system  (1).  This  verifies  the  stability  properties  or 
Integrity  property  of  regulators  of  the  special  class  of  regulators  described 
In  the  introduction. 


Characterization  of  Regulators  Defined  by  (2)  and  (3).  A  simple  method  of 
characterizing  such  regulators  is  to  relate  their  closed  loop  parameters 
to  the  closed  loop  parameters  of  optimal  linear  regulators.  This  method  can 
be  applied  to  second  order  systems  with  a  single  Input,  but  appears  to  be 
Intractable  for  general  systems.  For  the  simple  example  with 


a  >  0,  b  >  0 


08) 


Single  input  systems  also  demonstrate  that  the  set  of  controllers  defined  by 
(2)  and  (3)  is  a  conservative  estimate  of  the  set  of  controllers  with  integrity 
because  all  stabilizing  controllers,  for  single  Input  systems  possess  Integrity 
with  respect  to  Input'  outages. 


A 
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2 

the  closed  loop  characteristic  polynomial  Is  £  +  C^s  +  Cg.  Consider  the 
two  control  laws,  defined  by  (2)  and  (3),  and  uR  defined  by 

uR  *  -BTPrx  (19) 

with  PR  >  0  satlsylng  the  Rlccatl  equation 

PrA  +  ATPr  +  Qr  -  prbbtpr  *  0  (20) 

The  corresponding  coefficients  In  the  closed  loop  characteristic  polynomial 
are: 


CjL  *  b  +  (q22/2b)  +  (q11/2ab),  C2t  «  a  +  (q11/2a)  (21) 

C1R  "  tb2  +  q22R  ♦  2(  JaZ  *  q„R  C2R  -  /<?  ♦  <|1)R  (22) 

The  sets  of  possible  coefficients  for  these  two  control  laws  may  be  depicted 

In  the  (Cj,C2)  plane  as  shown  In  Figure  1.  The  set  for  uL  Is  a  subset  of  the 

set  for  uR.  The  lower  boundaries  of  these  two  sets  coincide  (the  line 

segment  C  *  a,  C,  >  b).  The  upper  boundary  for  the  uD  set  Is  the  segment 
c  '2  2  * 

of  the  parabola  C-j  «  b  +  2(C2-a)  with  Cj  >  b.  The  upper  boundary  for  the 

set  Is  the  line  segment,  b(Cj  -  b)  »  C2  -  a,  >  b  which  Is  tangent  to 


FIGURE  1.  Possible  Closed  Loop  Characteristic 
Polynomial  Coefficients 
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The  sets  of  possible  closed  loop  roots  may  also  be  determined  for  this  example. 
Figure  2  shows  these  sets  for  the  case  of  a  *  b  *  2.  The  negative  real 
axis  Is  contained  In  the  possible  sets  for  both  u^and  Ug.  The  remaining 
set  of  possible  roots  for  u^  is  a  small  subset  of  the  set  of  possible  roots  for 
Ug.  Although  these  results  appear  to  be  impossible  to  generalize,  let  us 
note an  Interpretation  of  the  sets  of  possible  roots  In  this  example  that  may 
be  generalized.  If  we  Introduce 


FIGURE  2.  Possible  Closed  Loop  Roots 

■  a 

a  positive  scalar  parameter  a  In  the  matrices  Q  and  QR,  say  Q  *  flQ  and 
Qr  ■  «Qr  and  consider  the  loci  of  closed  loop  roots  as  a  tends  toi  Infinity, 
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these  loci  tend  to  zeros  associated  with  Q  and  QR.  For  this  example  there  Is 
at  most  one  such  zero  and  It  must  be  real  and  negative.  In  the  optimal  case* 

a 

the  zero  Is  a  transmission  zero  associated  with  QR  and  its  magnitude  Is  arbitrary. 

In  the  case  where  u^  Is  obtained  via  (2)  and  (3)  the  magnitude  of  the  zero 
associated  with  Q  Is  restricted  to  be  less  than  b. 

In  this  example  the  root  locus  of  Interest  Is  the  locus  of  roots  of  the  polynomial 

p(s*G$  *  s2  +  2[b  +  (a/2ab)(qji  +  aqgg)]  +  a  +  (a/2a)  qjj  (23) 

As  «  tends  to  Infinity,  the  polynomial  may  be  factored  as 

p(s,a)  -  [s  ♦  bq„/(q„«q22)  ♦  0(<.'')][s  *  (<./2«t>)(q„«qM)  *  0(1)]  (24) 

so  that  one  root  tends  to  -bq^/(q^+aq22)  -  *  ^  an<^  r0°*  tends 

to  Infinity. 

Let  us  now  return  to  the  general  case  for  the  system  (1)  and  Introduce  the 
parameter,  a.  Into  the  control  law,  l.e. 

u  *  -aB^Px  (^®) 

with  P  given  by  (3V  Such  a  system  may  be  characterized  by  asymptotic  properties 
as  the  parameter,  o,  tends  to  infinity.  The  return  difference  matrix  Is 

T(s)  -  I  +oBTP(sI-A)“1B  (26) 

Algebraic  manipulation  leads  to  the  result  that 

TT(-s)  T(s)  «  I  +  oBT(-sI-AT)"1Q(sI-A)“1B 

♦  o2  BT  (-sI-AT)-1  PBBTP(*I-A)"1B  (27) 

Let  the  dimension  of  x  be  n  and  the  dimension  of  u  be  m.  For  simplicity, 
let  us  assume  that  BTPB  has  full  rank.  A  generalization  to  the  case  In  which 
the  rank  of  BTPB  Is  less  than  m  Is  of  Interest  but  Is  also  somewhat  more 
complicated.  If  we  let  s  *  ao  In  (27)  and  let  a  tend  to  Infinity,  we  obtain 
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TT(-oct)  T(wVI  -  ct'1  a"2  BTQB  -  <f2BTPBBTPB 
-►a-2  [o2I  -  BTPBBTPB] 


(28) 


We  also  note  that 

4>c(-s)4>c(s)  "  det  T(s^  (29) 

where  4»c(s)  denotes  the  closed  loop  characteristic  polynomial  and  $  (s) 
denotes  the  open  loop  characteristic  polynomial.  Since  the  closed  loop 
system  Is  stable  for  all  a  *  0,  we  can  deduce  from  (28)  and  (29)  that 
m  closed  loop  poles  tend  to  infinity  and  are  asymptotic  to 

s“  *  -a A1  (BTPBBTPB) »  1*1,2 . m  (30) 

m  -a  X^(B^PB),  1*1 , 2,..., m 

where  x(A)  denotes  an  eigenvalue  of  the  matrix  A.  The  remaining  n-m  closed  loop 
poles  approach  finite  values  which  are  the  left  half  plane  Images  of  the  zeros 
of  the  determinant  of  B^P(sI-A)~^B.  Denote  these  left  half  plane  finite  zeros 
by  s9,  1*1, 2,..., n-m.  The  eigenvectors  associated  with  the  finite  zeros  are 

1  T 

orthogonal  to  B'P  and  are  given  by 

X°  -  (sj  I -A)-1  Bv°  ,  1*1 ,2,..., n-m  (31) 

with  the  pj  determined  by 

BTP(s°  I-A)"1  B  yj  •  0.  1-1,2 . n-m  (32) 

The  eigenvectors  associated  with  the  asymptotically  Infinite  eigenvalues  are 
given  by 

X1  *  B  Vr  1,1  »2» ••••"» 


(33) 
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with  y*  determined  by 

Xj(BTPB)y“  *  -BTPB  y",  1*1 , 2,..., m  (34) 

These  results  are  similar  to  those  obtained  for  optimal  regulators.* 

For  the  optimal  regulator,  however,  the  asymptotic  eigenvalues  and  eigenvectors 
are  related  to  the  weighting  matrices  QR  and  RR.  Here,  unfortunately  the 
relation  Is  to  P  which  Is  In  turn  related  to  Q,  but  a  direct  relationship 
to  Q  Is  not  available. 


In  the  special  case  for  which  BT  *  [0  bJ]  with  Bj  an  mxm  nonsingular  matrix 
which  would  be  the  common  representation  for  systems  with  actuator  dynamics 
we  can  proceed  one  step  further.  In  this  case,  let  v"  *  Bjp“.  Then  from 
(33) 


» 

t 


s  ■[:,]  ’>  ■ 

and  from  (34) 

Xjvf  «  -B1(BTPB)  Bj“*  v"  *  -B^j  P4  vf.  1*1.2. (36) 

where  P^  Is  the  lower  right  mxm  block  of  P.  If  we  set  N  *  [v*»V2»***»vp  and  A  . 
*  dlag  (X^,  then  (36)  may  be  written  as 

NA  ■  -BjB^PjN  (37) 


0 

LV1 


,  1-1 ,2,...  ,m 


(35) 


Thus, 


P4  *  -(B1bJ)^HAIT1  (38) 

and  the  fact  that  P^  Is  symmetric  Implies  constraints  on  N  and  A,  which  may  be 
Interpreted  as  constraints  on  actuator  couplings  and  relative  bandwldths. 

•"Quadratic  Weights  for  Asymptotic  Regulator  Properties",  C.A.  Harvey,  G.  Stein, 
IEEE  Trans,  on  Auto.  Control,  Vol.  AC-23,  June  1978. 
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The  fact  that  the  finite  zero,  s®,  was  constrained  In  the  simple  example  Is 
a  property  that  Is  common  to  the  general  case.  Adding  2s®P  to  both 
sides  of  equation  (3)  and  rearranging  yields 


P(s®I-A)  +  (s°I-A)TP  -  Q  +  2s®P  (39) 

Multiplying  equation  (39)  on  the  right  by  V®  =  (s°l-A)  Bp^  and  on  the  left 
by  (v®)T  yields 


2(p°)TBTP  v®  =  (v®)T  (Q  +  2s®P)vJ 


But,  BTP  V® 


BTP(s®I-A)"1  By® 

(»°)T  Q 

2(v®)T  P  v® 


0  from  (32),  so  that 


(40) 


(41) 


and  It  follows  that 


(42) 


where  o(Q)  Is  the  largest  singular  value  of  Q  and  &(P)  Is  the  smallest  singular 
value  of  P.  Thus,  (42),  shows  that  the  magnitudes  of  the  finite  zeros  are  bounded, 
but  the  bound  Involves  P  and  Q.  Since  P  Is  a  function  of  A  and  Q,  the  bound  is 
an  Implicit  function  of  A  and  Q.  Unfortunately,  the  explicit  dependence  Is  not 
evident. 

An  alternate  characterization  of  controllers  defined  by  (2)  and  (3)  can  be  derived 
by  considering  the  optimal  controllers  for  (1)  with  respect  to  the  performance 
Index 


J  *  J  ($x^Qx  ♦  uTu)dt 
■'o 


(43) 
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with  0  a  small  positive  scalar  parameter.  In  this  case  the  optimal  controller 
can  be  represented  as 

u  a  “BT( j?  01  Pf)x  (44) 

1*6 


where 


( 


®1pi)  A  +  fi\)  +  BQ  * 


(Xb1  PJ  BBT(XB1P  ) 
1*0  1  1-0  p  V 


Equating  terns  of  like  powers  of  0  In  equation  (45)  yields 

poA  +  A\  “  p0  8  ®T  P0 

P,A  +  ATP^  +  Q  *  PqB  BTP1  +  P^  BTPq 


(45) 


(45a) 

(45b) 


etc. 


With  A  being  a  stable  matrix,  PQ  *  0,  so  that  the  right  hand  side  of  (45b) 
is  zero.  Thus,  the  equation  for  Pj  is  the  same  as  equation  (3),  and  controllers 
defined  by  (2)  and  (3)  may  be  viewed  ^multiples  of  first  order  approximation 
to  optimal  controllers,  i.e. 

u  *  -B^x  *  -BT”J  [j  JlB1Pi]x  (46) 

1*0 

Another  way  of  describing  this  characterization  Is  to  consider  the  controller 
given  by  (2)  and  (3)  with  Q  In  (3)  replaced  by  0Q.  Then  this  controller  is  the 
first  order  approximation  to  the  controller  given  by  (44)  as  0  tends  to  zero. 
This  Implies  that  the  closed  loop  root  loci  (parameterized  with  6)  associated 
with' these  controllers  are  tangent  at  0*  0  which  corresponds  to  the  open  loop 
roots. 


